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ABSTRACT 


Topography has long been measured by using stereo photographic 
image pairs collected from aircraft platforms. Space based stereo 
observations have not yet been fully exploited but interest in high 
resolution stereo coverage is increasing and the French SPOT satellite 
will begin to collect IS meter stereo pair Images of large areas of 
the eaith's surface in 198A. 

Conventional stereo coverage has been based on opposite side 
image combinations using the same sensor type. Same side offset 
views with identical sensor images have also been used. The Phoenix 
Corporation, under contract to the Jet Propulsion Laboratory, has 
been experimenting with quite a different technique for producing 
topographic information. The approach is based on same side/ same 
time viewing using a dissimilar combination of radar Imagery and 
photographic images. Common geographic areas viewed from similar 
space reference locations produce scene elevation displacements in 
opposite direction and proper use of this characteristic can yield 
the perspective information necessary for determination of base to 
height ratios. These base to height ratios can in turn be used to 
produce a topographic map. 

The experiment was performed using Synthetic Aperture Radar 
imagery from the Seasat Satellite and Return Beam Vidicon imagery 
obtained by the Landsat - 3 satellite. A tes’. area covering the 
Rarrisourg, Pennsylvania region was observed by these two systems 
in close time proximity and has the added value of having good co- 
planar ground control points in the common image area. 

The techniques developed for the scaling re-orientation and 
common registration of the two images are presented in this report 
along with the topographic determination data. The Harrisburg area 


Is well atspped and ics geographic and geologic classification has 
been studied extensively. Topographic determination based exclu* 
slvaly on the Images content Is compared to the nap Information 
which Is used as a performance calibration base. 


1.0 STUDY OBJECTIVES AMD SCOPE 


A major factor Influencing target spatial location Is the target 
altitude over rhe Image datum. Figure I, illustrates the direction 
and order of magnitude of the ground displacement of target position 
in the image because of altitude or topographic elevation. This 
displacement can be large, as indicated, with the direction and mag- 
nitude of the displacement a function of sensor type, l.e.. Active 
Sensor ■ AR (SEASAT-SAR) , Passive Sensor •• aR (Landsat RBV), look 

8 O 

angle and the range location of specific topographic areas. In Figure 

1, A Is the elevation of the image point T, and (a) Is the altitude 

of the satellite, S, with respect to the image datum plane. In a 

Radar-SAR- Image T will appear at and would have a ground range 

displacement of R . This displacement is opposite from that of a 
8 

passive optlcal-RBV-perspective, In which T would be Imaged at T^. 

It is interesting to note that If a target area were Imaged op- 
tically and with a radar, and properly correlated, all of the Informa- 
tion necessary to construct a three dimensional map of the area and to 
position ground targets at their correct map coordinates could be ob- 
tained. 

This was the objective of the Phoenix analysis using Images of 
the same geographic area produced by the SEASAT-Synthetic Aperture 
Radar and the Landsat-3 Return Beam Vldicon Imager. 

The Implementation of a Topographic Mapping Experiment required 
the development of a method for detection and correlation of common 
target points between the two images (geographic pair) to produce the 
final image products. 

The displacements shown in Figure 1 do not Include the effects of 
the earth's curved surface, which may be significant but are well known 
and were allowed for In study results. 

One potential source of error in the ultimate correlation of the 
image pair Is geometric misalignment due to the Inclination and re- 
sulting platform velocity vector difference between the SAR and RBV 
images. 
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Figure 1. Translation of elevated targets to displaced locations 
In Che Image plane. 



There ere also inege scale factor variaclons due to Che Inscan- 
taneoua Field of View differences between the SEASAI SAR and the RBV 
linages. 

In Che simplest case (i.e. , one in which there is no terrain relief 
and Che remotely sensed image returns are not a strong function of il- 
lumination aspect, time of date etc.) the effect of rotation and scale 
factor uncertainty is chat all pixels are not shifted the same amount 
so that the position offset is different throughout the scene. 

To the extent that this misalignment is knoira in advance, it can 
be determined and was compensated for by digitally resampling (i.e., 
rectifying) one reference Image to align it with the second image. 

In the studies performed by Phoenix this rectification was accomplished 
by using mapping equations derived by minimizing residual errors at 
analyst-selected control points. Bilinear polynomlnal equations can 
compensate for differences in scale, aspect ratio, rotation, skew and 
displacement between images assuming that there are no translations 
caused by heighting displacement. It is possible to select control 
points to minimize the heighting effects, (i.e., rivers, highway Junc- 
tions, airports, ...) and to calculate the control point heighting 
modulation on the alignment process when it occurs. The residual dif- 
ference at each control point and the rms value for all control points 
provide a measure of the registration accuracy. 

After registration, alignment, the true offset between the two 
Images is due to heighting translation and residual registration and 
scaling tolerances. The magnitude of this distortion is a function 
of : 

0^ SAR * Incidence Angle of the SAR Image to Particular Ground Points 
0^ RBV ■ Incidence Angle of the RBV Image to the Same Ground Points 
which are well known quantities. 

Using the relationship described in Figure 1 and calculations of 
Range based on orbit histories, a base to height ratio (relative) 


can be calculated and a topographic profile developed. This profile 
was developed and then compared to the actual map derived helghtlng 
translations which were used as performance baseline. 


2.0 CONCLUSIONS AND RECOWCNDATIONS 


Figure 2 suaBurizee the degree of eucceee eccelaed by ueing the 
image information to determine relief (helghting induced range mi* 
gration). The idealized or expected ahift waa calculated by uaing 
detailed map information and knowledge of aenaor(a) viewing incidence 
in the map area imaged. The incidence anglea were derived from pre- 
cise knowledge of satellite location and sensor pointing angles con- 
tained in attitude/orbit histories. Using this information the range 
(across the common imege area) translation versus helghting relation- 
ships were determined. The measured or Image(s) derived range of 
values were attained by first re-sampling the Images to a coiraor. 
orientation and pixel scaling base and registering both images to 
the same ground control point datum; then identifying point targets 
common to both images and measuring the image to image translation. 

Phoenix believes that the correlation between the expected and 
measured values is within tolerance considerations which Include: 

Variance - Expected Range of Values 

Incidence angle tolerance - minimal effect 
Translation (ground range) aspect - minimal effect 

Variance - Measured Range of Values 

Target recognition tolerance ± 2 pixels > 70' 

Registration tolerance ~ 85' 

Geometric Foreshortening————— - 100' 
in range (SAR Image) 

The variance in measured values is created by several uncertainties. 

One uncertainty is pinpointing the coonon targets in the RBV and SAR 
images. It is estimated that these uncertainties could be as much as 
2 pixels. However, target enhancements (contrast adjustments) were 
repeated many times to achieve the best possible edge on perimeter 
definition and it is felt in m.iuv cases that the resulting uncertain- 
ty was with 1 pixel and about 35 feet. The 35 feet pixel dimension 
is based on the re-sampled pixel dimensions. The uncertainty could 
have been in any of four directions in the planar image. 
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A second uncertainty are the co-registration tolerances. Regis- 
tration of each of the images was carefully done using a common set 
of control points. The results of the Interim or mid term report 
^ee Figure 3) revealed that comparison between expected values and 
measured values produced a fixed offset on bias between the results. 
Registration was redone using additional control points and making 
sure chat common control points were employed. The initial regis- 
tration used river area control points, as did the second registration, 
but due to the radiometric variability or diversity for common targets 
in Che two images all registration control points were not common to 
both images. The second registration corrected the systematic off- 
set but some tolerance or error in registration between the two Images 
still exists and is estimated to be about 2s pixels or 85 feet. As 
in the case of the target recognition uncertainty thr co-registration 
tolerance could occur in any direction in the planar image field. 


It was observed during the experiment that the SAR image was 
not map or ground range corrected. The Radar Instantaneous Field 
of Vlev( IFu^ in the range of cross-track dimension is: 

c 


IFOV 


2$ COSv 


c * The speed of light 

S « The bandwidth of the transmitted pulse (FM linear chirp) 

V => The target local grazing angle (radar vector intersection 
with the local tangent plane) 

The quantity is a constant value dependent on radar design para- 
meters and is referred to as the slant range IFOV. To produce a 
ground range corrected IFOV the must be used and the grazing 

angle V varies in range across the swath. The SAR image could be 
made Co line up with so-range features on the Harrisburg map but 
accurate correlation to map scale co-ordinates could not be achieved 
over Che extent of the range dimension. Ground control point regis- 
tration can partially correct for this condition but large numbers of 
control points throughout the image would be required to fully correct 
this anomaly. As a result an offset occurs between measured and ex- 
pected values and this offset has the effect of decreasing the mea- 
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sured values away from expected values with a clockwise aogular shift 
or in other words the tolerance is directional and not omnidirectional 
in the planar image field as were prior uncertainties. 

Slight variations in the SAR image were also noted in the along 
track or azimuth direction. This could be due to pulse records being 
dropped in the Doppler record prior to image signal processing or 
correlation. This would also produce a foreshortening that would 
also cause the unidirectional shift noted in the results. 

The RBV data matched map co-ordinates very well. 

The omnidirectional uncertainties combine to produce: 

Uncertainty - + 85^ 

- 110 ' 

This shift and the directional uncertainty of 110' account for 
some of the differences noted in Figure 2. Other factors such as the 
radiometric character of targets in each image, map uncertainties, 
attitude orbit error buildup may cause larger changes than expected 
and the degree of difficulty of precise registration make it a prime 
candidate for further work. 

In the authors opinion, the study results support the stereo 
potential of this unconventional approach and future efforts using 
this technique should be encouraged. The next opportunity to apply 
some of the techniques is with SIR-A experiment results using a sche- 
duled Landsat - 3 RBV image confluence. 


3.0 STUDY METHODS AND APPROACH 


3.1 DATA SELECTION 

The Important factors Influencing data selection were: 

1. To select images that were collected closely in time (days) 
to avoid seasonal or even weather (soil moisture, vegetation 
differences) related variances. 

2. To find an area that possessed a good control point population. 

3. To search for a moderate relief scene so that SAR translations 
not to be so severe as to interfere with the experiment. 

These criteria were applied after determining available candidates that 
merely complied with the more fundamental conditions of SAR and RBV 
image confluence in a cloud free area (a condition Imposed on the RBV 
not the SAR). The candidates that satisfied these more fundamental 
conditions were few in number, and the more selective criteria out- 
lined above narrowed the choice to the Harrisburg, Pennsylvania frame(s). 

The SAR image, during an ascending orbit, was collected by the 
SEASAT satellite on the 25th day of September, 1978 at 15:09:11 GMT. 

The RBV image was collected by the Landsat 3 satellite on the 29th 
day of September, 1978. Landsat 3 is in a sun synchronous orbit with 
a 9:30 a.m. equatlonal crossing time. This temporal adjacency complies 
with the study criteria to select images whose radiometric comparison 
would not be greatly masked by seasonal variances in target signatures. 
Study results indicate that this requirement was a necessary one as 
foliage, clear cut and other seasonal distinctions were very iiroortant 
to RBV and SAR specific point classif Icaitons. 

As shown in Figure 4 the SAR image was collected by an ascending 
SEASAT orbit with an inclination of 108° and an altitude of 799 kilo- 
meters. The Landsat - 3 RBV image was collected by a descending orbit 
with a 99° inclination and an altitude of 918 kilometers. The subimage 
used for the study is one scene - 99 km square that was acquired by 
one of the two RBV cameras. A complete RBV scene is a composite of 
two each 99 km x 99 km frames taken by two cameras. One camera col- 
lects data left of nadir and the other camera is pointing to the right 
of the nadir trace as the satellite follows its navigation path. 
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Figure 4 illustrates the relative orieotatlon and positioning of 
the RBV and SAR images. The image to image interactions that permit 
contour assessments are based on target/object range translations due 
to target elevation. As shovm in Figure 4 these translations differ 
in sense or direction and magnitude. The hypothesis used to conduct 
the study is sketched in the lower portion of Figure 4. The ground 
range translations in the two images are governed by the viewing per- 
spective. In the RBV image the translation in range AR is away from 
the sensor with magnitude: 

AR RBV - TAN9 x Object Height 

In the SAR image the translation in range AR is toward the sensor with 
magn) tude: 

AR-,_ - OBJECT HEIGHT 
TAN6 

in both cases 6 represents the target local angle of incidence. Just 
below the plan view of the RBV and SAR image frame overlays, the vector 
relationships and approximate translations for the study images are il- 
lustrated. 


3.2 RELEVANT SENSOR PERFORMANCE 

Table 1 summarizes some of the important sensor characteristics 
of the Landsat-3 RBV and SEASAT SAR. The RBV is a panchromatic camera 
with a ground IFOV of 24m. The RBV sensor produces an image using 
passive radiometric techniques while the SAR collects a doppler phase 
history that must be extensively processed with ground based signal 
processing equipment before an image is created. Although the manner 
in which these sensors create an image is very dissimilar some impor- 
tant characteristics, that are significant to the type of comparisons 
made during this experiment, have compatible values. The IFOV's are 
close together, ground coverage extent for some side viewing is near- 
ly equal and sensor signal to noise ratios are of the same general 
magnitude. 
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Figure 4 
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By concraac, cha radlomacrlc raspoasa to chana specific alamants 
in a acana is quits different for each of the sensors. A clear cut 
area can appear dark (low baekace'.:tar/intansity) in the SAR laaga 
and may appear bright (high reflected radiance) in the RBV image. 
Vegetation produces a relatively bright response in the SAR image and 
is darker in the RBV image. These and many more examples are contain- 
ed in the exhibit that serves as the major addendum to this report. 

Table 2 provides Important information about precise orbital 
and sensor pointing characteristics that must be used to accurately 
calculate the target local incidence angles in each of the scenes 
evaluated. Shown In Table 2 are the orbit and sensor pointing cha- 
racteristics for the SEASAT/SAR at the time that it observed the 
Harrisburg area. These quantities were used to calculate the range 
of range translation values expected across the SAR swath width and 
for specific locations where targets were singled out. 


3.3 STUDY METHOD 


The study method can be susnarized as follows: 

1. Several test areas common to both the RBV and SAR images were 
selected In the Harrisburg scene. Each of the test areas had 
interesting common features such as radio towers, beacons, 
power line crossings etc and the added criteria of possessing 
co-planar control points for common image registration. 

2. Each test area was re-sampled to accomplish a co-registration, 
a common RBV and SAR North realignment and a common pixel or 
picture unit scaling. After this was completed several con- 
trast enhancement or contrast stretch Intensity distributions 
were made for each of the candidate test areas to highlight 
the common features selected as point targets. 

3. Translation measurements were made and compared with the ex- 
pected translations using large sublmage frames (a few km's) 
to facilitate accurate measurements. 

Maps and image frames for each of the test areas are Included in 
the extensive exhibit that accompanies this report. Table 3 lists the 


14 - 


TABLE 2 SAR ATTITUDE ORBIT INFORMATION 
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1H.421 
18*421 
3H* 37* 
7V889V.02O 
1437*430 
-O.iOO 
-4V24.VV2 
9154*411 
179*144 
4729*841 
842*241 
840*907 
19*429 
24*949 
•211*197 
7184*099 
49.08S 


19.9I3 

19*247 

92*97 

4349*849 

0.1H7 

2489.094 

940*495 

189*741 

4729*472 

5*419 

9*741 

21*434 

24.977 

-711.712 

7183*871 

49*091 


^J49 

1*243 

280* 


743143* 

4374* 

0 * 

4499. 

5399* 

703. 

4779* 

851* 

870. 

23* 

24. 

- 212 * 

7183. 

49* 


19*913 

22.497 
0* 99*11 


217.442 

4729.337 

9.477 

9*809 

74.942 

24*987 

-212*894 

7183*417 

49.097 


D8N 

DSM 

PEGfHlHfSEC 

HETEkS 

KILOMETERS 

DEGREES 

KlLGNEIfRS 

HE7ERS/8EC* 

UNZTLES8 

HETERS/Sf.C 

KlLUNETERSfHlLlSCCnNnS 

K1LONETERS»NXL1SEC0ND8 

DEGREES 

KJLOHETERS 

HILIKCPINDS 

NETfRS/SEC* 

DEGREES 


SNATH 

%• 

39* 39* 94*71 

782* 

38* 

38*14 

DEGREES * HI NUTE8 » SECUNDS 

SNAIH 

2« 

jy* 41* 28*41 

282* 

94* 

31*79 

nESREES»HlNUTES»SECIINim 

KNAIH 

3*b 

39* 47* 0*43 

283* 

10* 

28*18 

DEGREES • MINUTES • SECflNbS 

SNAIH 

4 » 

39* 92. 30.14 

783* 

24* 

77*37 

DEGREEStHlNUTESpSECIlMDS 

SWATH 

9- 

39* 97. 97*70 

283* 

42* 

;!9*28 

DEGREES »HI NUTE8 »SECUNDS 


-3*983 

-3*281 

-7*444 

-8.798 


-3.829 

-3.749 

-8.191 

-8*734 


-J.289 

-4.423 

-8*949 

•8.497 


-3*018 

•9*381 

•8*741 

-8*839 


2.979 

«3*2S1 

DD 

4*48G 

-7*444 

D» 

8*731 

-8*798 

DB 

8.S99 

-9.149 

DB 


REN 1294 HIL 


Note: Incidence angle Is to near range side of 25 km segments of the 100 km 

swatch* 

^Important locations and incidence angle parameters 
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maps and ioaga frames chat comprise the exhibit. As indicated, six 
subimage test area subframes were finally selected after a search of 
detailed maps and the respective images for recognisable targets. Hie 
first map in the exhibit dlspLv.ys the entire Harrisburg area and the 
SAR- FRAME and RBV-FRAME orientation with respect to the map. 

The Phoenix corporation has developed a unique way of displaying 
images using dot-matrix printer techniques. The image pixel elements 
are portrayed using an array of dots (or absense of dots) to repro- 
duce grey levels. The Images created in this report use 17 grey 
levels. A special technique to prescribe (program) dot pattern 
printing and avoid systonatic noise or pattern effects has been de- 
veloped by Phoenix after extensive experimentation. The basic CCT 
tape assifflulation a'.d image production program also includes image 
rotation on re-ori.ritatlon (l.e.: North alignment) and pixel scaling 

(map matching or image to image common scalit.g) algorithms. Examples 
of these techniques and image registration are contained in the images 
that make up the exhibit of study results. 

3.4 DESCRIPTION OF TECHNIQUES USED IN EXPERLMENT ANALYSES 

The registration, re-orientation and scaling methods and procedures 
used to facilitate the study are briefly described in this section. Ad- 
ditional details can be supplied to interested readers upon request. 

The procedures outlined herein were implemented on the Phoenix Corpo- 
rations Prime 550 computer system. 

Figure 4 provides an illustration of how the RBV and SAR Images 
cover the Harrisburg area. The SAR nadir trace is a considerable dis- 
tance from the image (~ 300 km to the image centerline) while the RBV 
nadir trace is in the center of the frame. The inclination differences 
account for the general frame orientation shown. 

Using spherical trigonometry the offset angle (discounting attitude 
differences) can be calculated as follows: 
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LOCAL MERIDIAN 



6 ANGLE BETWEEN 
ASCENDING OR DESCENDING 
PASSES OF A SATELLITES 
ORBIT WITH RESPECT TO 
THE LOCAL MERIDIAN. 


A, 

- TM - 

>> 2-. 
cos 9 ) 

z 

*^080 

2 J 


- cos 1 


SEASAT 

LANDSAT - 3 

N 

3 

18 

R 

43 

252 

1 

108° 

1 ,0 


1 ■ Orbital Inclination 

6 " Local latitude 

N “ Number of days before 
orbit begins a repeat 
cycle. 

R ■ Number of orbits before 
a repeat cycle begins. 

Harrisburg, Pennsylvania 

Latitude ■ 40° 17' 


* Figure 3. Orbital crossing angle 
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SEASAT CALCULATION 


A, - TAN-^ 
z 


-.309 « .0406 



26.4® 


LAtn>SAT CALCULATION 



la ch« accuaI rA-orlAnCAClon progrAms ground control point locations in 
AAch inage were used to implement a common image orientation. 

Registration 

The area around Harrisburg^ Pennsylvania, visible in both *:he SAR and 
RBV scenes is dominated by the Susquehanna River along with several small- 
er tributaries and streams. Within the river there are numerous islands 
and other features clearly distinguishable on the SAR image. Immediately 
adjacent to the river are several major highways and interchanges visible 
on the RBV image. These features provided a source 'it control points on 
a level surface in an area that is otherwise dominateu by variations in 
relief . 


Preliminary inspection of the two images revealed that it would be 
difficult to find one-to-one correspondences between features or indi- 
vidual pixels between the two images. However, each image showed r^tadlly 
recognizable features that could be identified on topographic maps. A 
decision was made to register each image separately to a geodetic co- 
ordinate system using the clearly seen features in each image and the 
corresponding map location. Special care was taken to insure that each 
point selected for control was on the same level surface. Control point 
locations were shown in the mid term report. Figure 4 indicates that 
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expected translation to elevation ratios are nominally 0.1 for the RBV 
scene and 2.0 -f 2.8 (far to near range) for the SAR scene dependent 
on target location. Hence, nominal elevation differences between con- 
trol point surfaces a few meters apart will result In translation of only 
several meters; far less than the size of a pixel. 

Each of the areas selected for analysis referenced In table 3 were 
separately registered and resampled using the same control points. The 
actual registration was done in the following way. Each study area was 
defined by its latitude and longitude perimeter extent, designated 

LATMIN 

LATMAX 

LONMIN 

LONMAX 

Within each individual area the latitude and longitude wore recomputed with 
respect to LATMIN AND LONMIN where 

LAT' » LAT - LATMIN 
LON’ - LON - LONMIN 

For each control point, a LAT and LON value was obtained from the topogra- 
phic map and a corresponding row and column location I, J was obtained 
from the respective SAR and RBV images. These matrix and geodetic loca- 
tions were used in a least squares solution to find the coefficients to 
the following equations. 

I - Cl * LAT' + C2 * LON' + C3 
J - C4 * LAT' + C5 * LON' + C6 

This form allows the image data to be rotated to North, allows for a 
translation between the new geodetic origin LATMIN, LONMIN and the row, 
column origin of the image data and also allows different scale changes 
along the orthogonal axes. It does not allow more sophisticated rubber- 
sheet capabilities such as skewing or warping of the image. It was felt 
that most of the gross effects would be removed during the original satel- 
lite data processing and in addition each individual area was so small 
that these effects would not amount to much if they were present, its 
indicated in the Section (2.0), the SAR range data was not corrected on 
a pixel by pixel basis for precise slant range to real range geometric 
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fidelity. This was discovered at the conclusion of the study when compar- 
ing measured and expected results. Better registration accuracy would be 
achieved by resampling SAR CCT's and correcting for this deviation. 

The success of our registration algorithm used In this study Is 
estimated from the rms deviation obtained In the following way 

. - 1 )^ 

°1 N 

Where 1' Is the recomputed value of the control points row location and 
N Is the number of control points used. A similar expression was used 
for Oj. The values of and are summarized for each Individual 
area In Table 4. It Is seen that the registration accuracy estimated 
this way was typically on the order of one or two pixels. Given that 
we would be looking at elevation differences of about 1,000 feet and 
corresponding translations of about 2,500 feet or about 70 resampled 
pixels this registration accuracy was felt to be acceptable. A re- 
sampled pixel Is 11.12 meters or 36.5 feet. The Table 4 representa- 
tion Is In original RBV SAR IFOV/Plxel units. Comparable 1,000 feet 
translations are about 30 pixels. In the conclusions section a con- 
servative registration tolerance of 85 feet was assumed. 

Resampling 

In order to measure displacements due to elevation using the SAR 
and RBV Images, It Is necessary Chat both Images be Identical In scale 
and orientation. The previous section explained the basis by which both 
Images are oriented to north. The actual procedure during which this 
is carried out, while at the same time a scale size is specified to pro- 
duce a normalized pixel, Is called resampling. The resampling Is a map- 
ping process by which a grid containing the original data whose row 
and column orientation and grid spacing Is characteristic of the Indi- 
vidual sensor Is recomputed to some other grid orientation and spacing. 

The original data sets were resampled for each of the six test areas. 
The resampling limits were defined by LAITIIN, LATMAX, LONMIN, LONMAX va- 
lues. The resampling began at the location LATMIN, LONMIN proceeded at 


- 22 - 


TABLE 4 


Summary of Registration Accuracies 


Area 

Si 

\R 

RBV 





1 

0.9959 

1.7337 

1.5112 

0.8074 

2 

0. 9958 

1.7336 

1.5111 

0.8077 

3 

0.9951 

1.7340 

1.5125 

0.8077 

4 

0.9960 

1.7339 

1.5116 

0.8071 

5 

0.9959 ! 

1.7336 

1.5113 

0.8078 

6 i 

0.9956 ' 

1.7338 

1.5121 

1 

0.8078 
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constant latitude from LONMIN to LONMAX. The resampling latitude was 
Chen Increased by A LAX and the resampling repeated. The sampling In- 
terval In longitude, A LON, was defined by 
A LON + A LAT/COS (LAT) 

This ensured Chat equal areas were sampled In the east-west and north- 
south directions. 

For each LAT, LON position In the resampling grid a corresponding I, 
J location was found In the Image data using the same coefficients solved 
for In Che the registration algorithm. From this I, J value the four 
image pixels surrounding the LAT, LON position were identified and used 
to linearly Interpolate an amplitude at the LAT, LON position. This 
procedure Is Illustrated In Figure 6. 

The same latitude and longitude limits were used for each study 
area for both the SAR and RBV images. The latitude sampling interval 
was 0.0001 degree for both Images. This Interval yields a linear 
sampling distance that is slightly smaller Chan either the SAR or RBV 
ground sampling distance. This resampling and registration allowed 
the translation of features due to elevation changes Co be observed 
and measured by a direct comparison of the two Images. 

Translation 

Changes in elevation are perceived differently by active and passive 
sensors. This Is Illustrated In Figures 1 and 4 where a vertical height, 
h, Is viewed by sensor a large dlstarce away at an angle 6 with respect 
CO the local vertical. 

Since most elevations will be very much smaller than the distance 
to Che satellite borne sensor and smaller chan the radius of the earth, 
the representation in Figures 1 and 4, showing a locally flat earth, is 
perfectly satisfactory. However, the assumption of a flat earth cannot 
be made with regard to the angle 6. Due to earth curvature, this angle 
can change within a field of view and in general Is not equal Co Che 
antenna or optics viewing half angles. Figure 7 shows how 6 may be 
computed. 

Consider a spherical earth of radius r^, a reasonable assumption 
over the distances being considered. Consider also a sensor at alcl- 


- 24 - 



Figure 6. The relation between a point LAT, LON in the resampling 
grid and the four surrounding points In the image grid 
used for Interpolating a value at LAI, LON. 





tude H whose subsacelllte location Is at a latitude and longitude given 
by 6^, X^4 This sensor views a height h at location 6, X, a distance d 
away. Constructing the line z we see that 

6 “ sin ^ {£) 
d 

We also note that 

2 - (r + H) sin 4' 
e 

d- + + 

where h « r _ 

and 4* ” Cos ^ I sin 0 sin 0 + cos 6 

L ° 

For a given sensor location, we are now able to compute the expected 
translations due to elevation as a function of object location. 

To note the magnitude of this effect consider an active sensor with 
an angle, 0', of 20° (depression angle > 70°). The flat earth approxi- 
mation would set 0 * 20° and compute a translation for a 100 m elevation 
of 27A.7 m. Such a sensor at an altitude of 800 km would be viewing 
objects displaced an angular distance 4' » 2.5°. Using r ■ 6371 km 
we find for this case 0 * 21.5°. This gives a translation of 253.9 m 
for the same 100 m elevation. Hence, for a system like SEASAT, a flat 
earth approximation would yield translations In error by about 8%. 

Inclination angles were calculated using target local geometries 
and the relationships described above. These relationships are shown 
In Figure 8. 

3.5 MEASURED RESULTS 

Table 5 lists the measured results obtained from each of the SAR 
and RBV target subframes contained In the exhibit. The x and y dis- 
tances to each target were accurately measured from the lower left 
corner of each subframe which serves as a common registration datum 
point. The location of the target In the SAR frame was determined and 
the target/local Inclination angle calculated which specifies the trans- 
lation value (2.0 to 2.8) and this was used to produce a AR AE relation- 
ship In Figure 2. Figure 2 reproduced as parr of this section. In- 


- 2 r (r + H) cos 4' 
e e 

cos 0 cos 
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Figure 8- SAR KfiV &E/AR Relationships 
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Scale 652’ per inch longitude 


dividual target locations are shown in the reproduced version of the 
figure. Accounting for the tolerance variations and SAR range distor- 
tions described in Section 2.0 the comparison is technique validating. 

A similar experiment should be conducted using the upcoming SIR- A 
information. One or more of the principle investigation areas planned 
for SIR-A should be designated for the repeat analysis. 
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